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Abstract. STAR has measured a variety of strange particle species in p + p 
collisions at ^/s = 200 GeV. These high statistics data are ideal for comparing 
to existing leading- and next-to-leading order perturbative QCD (pQCD) mod- 
els. Leading-order (LO) models such as PYTHIA need to be tuned to describe 
identified strange particle data from STAR. We show that tuned PYTHIA 
can also describe the pt spectra of strange resonances. More rigorous Next-to- 
Leading order pQCD calculations using parameterized fragmentation functions 
for quarks and gluons will also be compared to STAR data. The OPAL ex- 
periment has recently released e'^ + e~ data from light quark flavor tagged 
analyses allowing for the first time to make precise parameterizations of light 
flavor separated fragmentation function. We show that our A data put a more 
stringent constraint on the gluon fragmentation function than e+ + data. 
Furthermore we show that pQCD fails to describe the observed enhancement 
of baryon-to-meson ratio at intermediate pt (2-6 GeV/c), which may be a first 
indication of other, non-perturbative mechanisms at play in p + p collisions at 
that momentum. 

Keywords: p + p collisions, pQCD, PYTHIA, fragmentation functions 
PACS: 25.75.Dw, 25.40.Ep, 24.10.Lx 



1. Introduction 

Perturbative QCD has proven to be successful in describing inclusive hadron produc- 
tion in elementary collisions. Within the theory's range of applicability, calculations 
at next-to-leading order (NLO) have produced accurate predictions for transverse 
momentum spectra of inclusive hadrons at different energy scales El- With the 
new high statistics proton-proton data at y/l — 200 GeV collected by STAR, we can 
now extend the study to identified strange hadrons as well as strange resonances. 
The perturbative QCD calculation applies the factorization ansatz to calculate 
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hadron production and relics on three ingredients. The non-perturbative parton 
distribution functions (PDF) are obtained by parameterizations of deep inelastic 
scattering data. They describe quantitatively how the partons share momentum 
within a nucleus. The second part, which is perturbatively calculable, consists 
of the parton cross-section amplitude evaluated to LO or NLO using Feynman 
diagrams. The third part consists of the non-perturbative Fragmentation functions 
(FF) obtained from e+ + collider data using quark-tagging algorithms. These 
parameterized functions are sufficiently well known for fragmenting light quarks, 
but less well known for fragmenting gluons and heavy quarks. Recently, Kniehl, 
Kramer and Potter (KKP) have shown that FF are universal between + e~ and 
p + p collisions [ '51 . 

In the following section we compare our p-fp data to PYTHIA, the most com- 
monly used leading-order Monte Carlo event generator for elementary collision. We 
then move on to compare our data with more sophisticated NLO calculations. This 
will lead to a discussion of the difference between quark and gluon jets in p + p. In 
the final section we discuss the different contribution of these two jet types to the 
production of baryons and mesons. 

2. Data Analysis 

The present data were reconstructed using the STAR detector system which is 
described in more detail elsewhere [ 0] . The main detector used in this analysis 
is the Time Projection Chamber (TPC) covering the full acceptance in azimuth 
and a large pseudo-rapidity coverage (| ?7 |< 1.5). A total of 14 million non-singly 
diffractive (NSD) events were triggered with the STAR beam-beam counters (BBC) 
requiring two coincident charged tracks at forward rapidity. Due to the particulary 
low track multiplicity environment in p-l-p collisions only 76% of primary vertices are 
found correctly; from the remainder, 14% are lost and 10% are badly reconstructed 
as a MC-study showed. Of all triggered events, 7 million events passed the selection 
criteria requiring a valid primary vertex within 50cm along the beam-line from 
the center of the TPC detector. The strange particles were identified from their 
weak decay to charged daughter particles. The following decay channels and the 
corresponding anti-particles were analyzed: Kg 7r+ -f 7r~ (b.r. 68.6%), A —>■ 
p + 7r-(b.r. 63.9%) ^ A + 7r-(b.r. 99.9%). Particle identification of the 
daughters was achieved by requiring the dE/dx to fall within the 3cr-bands of the 
theoretical Bethe-Bloch parameterizations. Further background in the invariant 
mass was removed by applying topological cuts to the decay geometry. Corrections 
for acceptance and particle reconstruction efficiency were obtained by a Monte- 
Carlo based method of embedding simulated particle decays into real events and 
comparing the number of simulated and reconstructed particles in each pT-hin. 
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3. Comparison to PYTHIA 

One of the most widely used models for simulating elementary collisions is PYTHIA 
It is a parton-shower based event generator that includes leading order parton 
processes and parton fragmentation based on the Lund Model. The parton dis- 
tributions of the initial state protons can be chosen from an array of PDFs (here 
we use CTEQ5M). The model is being actively used and the authors have recently 
released a version with completely overhauled multiple scattering and shower algo- 
rithms (version 6.3) [ 6 . The PYTHIA version used in this paper is 6.317. 

The string fragmentation based on the Lund Model requires only two parame- 
ters to define the shape of the fragmentation function and is universal for all light 
quark flavors. Baryons are produced from di-quarks and their probability is sup- 
pressed with respect to qq pair production. Next-to-leading order processes can 
be "simulated" in PYTHIA by tuning the K-factor (MSTP(33)) or by increasing 
the parton shower activity. This will enhance the relative probability of hard pro- 
cesses of type quark-gluon and thus mock-up the contributions from higher order 
processes. In figure 1 we first compare the measured STAR spectra for identified 
pions and protons to a simulation from PYTHIA. The pions agree very well with 
the default parameters whereas the protons seem to lie in between the default and 
the tuned K-factor calculation. In figure 2 (upper row) we compare PYTHIA cal- 
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Fig. 1. Identified 7r+ and protons in minimum-bias p + p collisions at \fs= 200 
GeV compared to predictions from PYTHIA v6.3 with and without K-factor. Data 
from [Q 



culations for strange mesons and baryons to the measured STAR data. Whereas 
the default parameters agree quite well for the Kg, they clearly underestimate the 
yields at intermediate pT for the A and By increasing the K-factor to 3 we 
achieve a reasonable agreement with the data. In figure |21 (lower row) we compare 
PYTHIA to the strange resonances K* , (p and S*. Again, only when applying a 
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Fig. 2. (Top) Minimum-bias pt spectra for Kg, A and S at (| y |< 0.5) from p + p 
at -^i = 200GeV. (Bottom) K* and (f>, and E* spectra at mid-rapidity. In the 
left panel blue squares are K'^* and black symbols are K'^* . Resonance data from 

[0 M nm 



In summary, PYTHIA is capable of describing px spectra for a variety of par- 
ticles from p + p collisions at RHIC energies. However, we have presented evidence 
that a tune of the LO K-factor is necessary in particular for strange baryons and 
resonances. Of course, we have not explored all possibilities of parameter "tunes" 
and there may be other, equivalent ways of reproducing the data. 

What are the possible reasons for this discrepancy? The "naive" reason, sup- 
ported by the K-factor tune, is that higher order contribution may be significant. 
However it is troubling that the pions in figure ^ do not seem to require this tune, 
thus introducing a rather "unnatural" particle species dependance. Nevertheless, a 
similar study of K-factors for non-identified hadrons found that at ^/s — 200GeV 
a value of 3 was needed [ II . Another, maybe more plausible explanation may be 
related to fragmentation functions in PYTHIA. It could be that some flavor depen- 
dant refinements to the Lund symmetric string fragmentation are necessary. This 
will be discussed in the next section with the use of NLO calculation and a set of 
parameterized FF. 
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4. Comparison to next-to-leading order pQCD 

A next-to-leading order calculation can help solve both these problems by including 
higher order parton processes and rigorously parameterized fragmentation func- 
tions. Fragmentation functions for separated quark flavors have been notoriously 
difficult to obtain due to the lack of sufhciently precise collider data. However, re- 
cently OPAL has published flavor tagged data which allowed theorists to compute 
better fragmentation fimctions [ I12j . 




Pj [GeV/c] Pj [GeV/c] 



Fig. 3. pt spectra for Kg (left), A (right) at midrapity (| y |< 0.5) bom p + p at 
^/s = 200GeV compared to two different NLO calculations. Dashed lines indicate 
the scale uncertainty of the NLO calculation, ie. fi — 0.5pT (lower), /i = 2pT 
(upper). 

In figure 13 we compare two different NLO calculations to our Kg and A data 
. The first one (black lines) uses older FF by Kniehl et al. (KKP) and Vogelsang 
et al (WV) [^5. The second one (red hues) was done by Albino et al. (AKK) 
using more recent FF based on light flavor tagged OPAL data Clearly, these 

newer parameterizations improve the description of our A data greatly. However, 
in order to achieve this agreement they fix the initial gluon to A fragmentation 
function {D^) to that of the proton, then estimate that a additional scaling factor 
of 3 is necessary to achieve agreement with STAR data. This modified FF for 
however also works well in describing the p + p SPS data at y/s = 630GeV. So it 
appears that the STAR data is a better constraint for the high z part of the gluon 
fragmentation function than the OPAL e+ + data. Similar conclusions with 
respect to the important role of p + p collisions have been drawn elsewhere [ llHj . 
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5. Hadrons from quark vs gluon-jets 

In order to investigate the fragmentation of gluons further, we looked at other 
observables that may be sensitive to this. For this study we used the PYTHIA 
Monte Carlo generator to differentiate events with gluons jets vs. quarks jets in the 
final state. We define a "Gluon-jet" event as one where the final partons are g-g or g- 
q and a "Quark-jet" event one where the final partons are q-q. The overall weighting 
of these events in the total p + p cross-section is dominated in favor of gluon-jets in 
PYTHIA. We then compared PYTHIA to the arbitrarily scaled ttit spectra from 
data. The scaling was adjusted such that all spectra would overlap in the 0-1.5 GeV 
region. 
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(a) Scaled data 
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Fig. 4. (Left)Arbitrarily scaled uit spectra for baryons and mesons from p + p 
collisions at -^/s = 200GeV. (center) Scaled ttit spectra for Gluon-jet events from 
PYTHIA. (right) Scaled tot spectra for Quark-jet events from PYTHIA. 

It is interesting to observe that gluon jets will fragment very differently into 
baryons and mesons than quark jets. For gluon jets, there is a clear shape differ- 
ence between baryons and mesons at niT ~ l.SGeV, consistent with the di-quark 
suppression in the string fragmentation picture. For quark jets, the shape differ- 
ence is modified by an additional dependency on mass of the produced particle. 
When comparing the result from PYTHIA with our data, this picture indicates the 
dominance of gluon jets in p -|- p at RHIC energies. 



6. Baryon production in pQCD 

In string models (i.e. PYTHIA) baryon production is understood via the production 
of di-quarks pairs from string-breaking and their recombination with other quarks. 
This process is suppressed with respect to q — q pairs from string- breaking resulting 
in systematically lower baryon yields than mesons. 
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(a) STAR data vs PYTHIA 



(b) STAR vs UAl 



Fig. 5. Ratio of A over Kg vs pt at midrapity from p + p at y/s = 200GeV and 
y/s = 630 GeV compared to predictions from PYTHIA. UAl data from [I17|. 



Recent heavy ion data from STAR show a large enhancement of the baryon to 
meson ratios at intermediate pr, which is associated with parton coalescence and 
recombination models [El- Figure [51(a) shows the predictions from PYTHIA for 
the A /Kg ratio vs px- We have separated the results for gluon and quark jet events 
to show how these contribute to the ratio differently. The overall PYTHIA result 
lies in between, closer to the red line. We observe that the prediction underestimates 
our data by at least a factor of 2 in the pr range from 0.8-3 GeV/c. 

In figure |31(b) we show that the disagreement is not specific to our energy scale 
but also exists at Tevatron energies. At ^/s = 630 GeV, the difference between 
PYTHIA and data is about a factor of 3 and the enhancement of A /Kg is twice as 
large as in STAR. This may be an indication that the effects observed in this ratio 
in heavy-ion data are present in some form in p + p data and therefore this should 
be noted before attributing heavy-ion phenomena to explain this effect. 



7. Summary 

We have shown that the theoretical description of identified strange particles in 
p + p and p + p collisions is still not fully understood. This is especially important 
since these models are now extensively used to predict observables for the LHC- 
era, and therefore one should be aware of their limitations. Phenomenological LO 
models can be tuned to describe the data but are unable to describe the baryon 
enhancement at intermediate px- NLO calculation have greatly improved with light 
flavor tagged fragmentation functions. However the high-z range of the gluon FF 
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previously extracted from + e~ data seems inconsistent with p + p and p + p 
data, indicating that RHIC data could be valuable in constraining the gluon FF. 
Arbitrarily scaled rriT spectra for strange particles exhibit mx scaling and confirm 
the dominance of gluon jets in p + p and therefore the importance of understanding 
gluon fragmentation. Finally, the baryon to meson ratio at intermediate px is 
about a factor 2 larger than predicted by pQCD. This difference is even larger at 
^/s — 630 GeV in p + p collisions. This is an indication that the baryon/meson 
effects previously observed in heavy ion collisions are present in some form in p + p 
data, and that the associated physics phenomena therefore need to be explained 
without requiring the presence of a quark-gluon plasma. 
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